function [y, T] = dynamic_3(y, x, params, steady_state, sparse_rowval, sparse_colval, sparse_colptr, T)
  y(19)=T(1)*(1-params(4))*y(15)^(-params(4));
  y(18)=T(4)*params(4)*y(3)^(params(4)-1)-params(3);
  y(16)=y(11)/y(15);
end
